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Abstract 

We provide an alternative view of the efficient classical simulatibility of fermionic linear optics in 
terms of Slater determinants. We investigate the generic effects of two-mode measurements on the Slater 
number of fermionic states. We argue that most such measurements are not capable (in conjunction with 
fermion linear optics) of an efficient exact implementation of universal quantum computation. Our 
arguments do not apply to the two-mode parity measurement, for which exact quantum computation 
S : becomes possible, see Q. 

Ph. Dedication to Asher Peres 
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It is a pleasure to contribute to the Festschrift for Asher Peres' 70th birthday. To characterize Asher's 

3 : 

achievements in physics and quantum information theory, we would like to quote from a novel by the 



X 



Austrian writer Robert Musil (1880-1942) The man without qualities written in 1930 0. The protagonist 
of the novel, Ulrich, is full of praise about science: 

But one thing, on the other hand, could safely be said about Ulrich: he loved mathematics 
because of the kind of people who could not endure it. He was in love with science not so much 
on scientific as on human grounds. He saw that in all the problems that come within its orbit, 
science thinks differently from the laity. If we translate "scientific outlook" into "view of life," 
"hypothesis" into "attempt," and "truth" into "action," then there would be no notable scientist 
or mathematician whose life's work, in courage and revolutionary impact, did not far outmatch 
the greatest deeds of history. The man has not yet been born who could say to his followers: 
"You may steal, kill, fornicate - our teaching is so strong that it will transform the cesspool of 
your sins into clear, sparkling mountain streams." But in science it happens every few years 
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that something till then held to be in error suddenly revolutionizes the field, or that some dim 
and disdained idea becomes the ruler of a new realm of thought. Such events are not merely 
upheavals but lead us upward like a Jacob's ladder. The life of science is as strong and carefree 
and glorious as a fairy tale. And Ulrich felt: People simply don't realize it, they have no idea 
how much thinking can be done already; if they could be taught to think a new way, they would 
change their lives. 

Asher Peres is one of those scientists whose work exemplifies the force of logical thinking and of in- 
dependent and unorthodox investigation into the nature of physics, as demonstrated by his wide ranging 
publications from relativity theory to quantum mechanics. We are grateful for his "revolutionary and coura- 
geous acts" in quantum information theory which have been an inspiration for our own work. We also hope 
that Asher will appreciate this small paper whose subject is at the boundary of physics and information, a 
boundary that Asher has enjoyed crossing during his productive scientific career. 

1 Fermionic linear optics and single-mode measurements 

A short introduction to second quantization 1 will serve to set our notation. Suppose that the Hilbert space 
of a single electron has dimension D, and is spanned by a standard basis | i), 1 < i < D. The label i may 
indicate both spin and space degrees of freedom. In the language of second quantization these same basis 
vectors are indicated as a\ \ 0). | 0) is the basis vector of the Hilbert space corresponding to no electrons 
(the vacuum state). a\ is the creation operator for an electron in mode (or orbital) i. Without the dagger it 
is a destruction operator. Creation operators are taken to anticommute 2 which enforces the Pauli principle 
since (a\) 2 = 0, i.e., two electrons cannot be put in the same state. 

We are interested in fermionic states that contain not just one electron but N < D electrons. In second 
quantized language, an example of such a state is 

nj^ajio). (1) 

This state is special not because it places electrons in the standard-basis orbitals, but because it can be 
'Apparently, the second-quantized analysis of fermions developed over the course of a series of papers, principally by P. Jordan; 
the most notable in this series is P. Jordan and E. P. Wigner, "About the Pauli exclusion principle", Z. Physik 47, 631 (1928). It had 
assumed essentially its modern form by the time of V. A. Fock, "Configuration space and second quantization", Z. Physik 75, 622 
(1932). It is reviewed in innumerable modern texts, e.g., A. L. Fetter and J. D. Walecka, Quantum theory of many-particle systems, 

(McGraw-Hill, New York, 1971), Chap. 1. 

2 {a\,a t j } = {a,i, aj} = and {ai,aj} = 6ijl. 
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written as a single term. In the old language of electron physics, Eq. Q is an example of a single Slater 
determinant. A general iV-electron state is a superposition of such Slater determinants, in which the product 
over orbitals can run over any set of iV orthogonal single-particle basis states. 

Recent work on quantum information processing with fermions has focussed on this Slater-determinant 
characterization of electron states. For example, a proposed measure of entanglement of fermi systems 
0IU is the "Slater number", the number of terms in the expression for the wave function involving the least 
number of Slater determinants. 

It has been known since the very earliest work in computational physics that the simulation of physical 
properties of electron systems becomes tractable if the states can be approximated by single Slater deter- 
minants (this is an essential feature of the Hartree-Fock approximation). This approximation is exact if the 
Coulomb interaction between electrons can be ignored. In fact this is a rather drastic approximation, and 
much of the art of atomic (and solid state) modelling has consisted of finding well chosen "mean fields", 
atomic potentials that mimic as well as possible the average effect of the interaction of many electrons in the 
atom. This endeavor has been rather successful, and has provided a basis for the approximate computation 
of many quantum properties in atomic, molecular, and solid-state physics, as well as in chemistry. 

Recent work by the authors [5] (see also Ref. [6|) has shown that the restricted quantum computational 
process of fermionic linear optics can be simulated efficiently on a classical computer. Fermionic linear 
optics on a set of non-interacting electrons are operations such as beam splitters, phase-shifters (delay 
lines), von Neumann measurements of the electron state, with the choice of quantum operations potentially 
based on prior measurement results. 

It was not explicitly shown in |5 1 that the Slater number remains one under these operations, and we will 
show it here and argue that, perhaps not surprisingly, it provides the basis for the classical simulatibility. 
First, the allowed Hamiltonian evolutions in this model are in the class of "one-body interactions"; that is, 
they arise from forces between the electrons and the controlling apparatus, and not between different elec- 
trons. Such Hamiltonians Hi have non-zero matrix elements only between iV-particle Slater determinant 
states $i and $2 with the same N, i.e., the Hamiltonian is number conserving (although an extension to 
fermion-parity conserving Hamiltonians is possible, and has been worked out in OHO); N—lof the orbitals 
in $1 and $2 should be identical, and just one may be different in the two states. For example, generally 

<$! m $ 2 ) ^ if I <&!> = Hj^aJl 0), I $ 2 > = (n^aj)^! 0) (2) 

but 

<<&! m $ 3 ) = if I $1) = n^ajl 0), I $3) = (n^ajK^aj^l 0) (3) 
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In this last example, the matrix element would be nonzero if two body terms in the Hamiltonian (electron- 
electron interactions) were included. In general, such a non-interacting Hamiltonian can be written as 



D 



Hit) = M.' ;»:<',• (4) 

Equation © introduces the Hermitian time-dependent matrix b(t). We will use the standard result, reviewed 
in [5 1, that the action of the time evolution operator, 



U(r)=Texp(-i dtH(t)) (5) 
Jo 

(T is the time-ordering operator), on a creation operator is 

U{T)a\u\r) =Y^V im (r)al = oJ(r). (6) 
Here the unitary matrix V is given by 

V(r) =Texp(-i f dth(t)). (7) 
J o 

The notation introduced in the last part of Eq. © (the r dependence) indicates that the resulting operator is 
just the creation operator for an electron in the time-evolved orbital 

\i(r)) =J2 V im(r)\m). (8) 

Under U (r), then, the initial state Eq. (Q evolves to (using U\ 0) = | 0)), 

n^aJ^lO), (9) 

i.e., still a single Slater determinant in a rotated basis. 

We now turn to the other computational step considered by |5 1, projective measurement of the occupa- 
tion of a single orbital (call it | k)). The projector corresponding to the state being occupied is 

Pi = a[a K , (10) 

and for the unoccupied outcome, the projector is 

P = 1 - a) K a K = a K a\. (11) 

What is noteworthy is that both projectors consist of a single product of annihilation and creation operators, 
which would not be the case for bosons. 
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Let us now show that the state after measurement, under all circumstances, continues to be a single 
Slater determinant. That is, we show that 

PoCD^xoJCrJJIO) (12) 

and 

P 1 (Ug 1 a\(T))\0) (13) 

are single Slater determinants, of a very simple form. The technical steps are described in an Appendix 
of E), we give a simple version of them here for completeness. We first write the orbital | k) as a linear 
combination of an orbital in the space spanned by the set {| i(r)}}, and an orbital not in that span: 

\ K ) =a\ in) + 0\ out) (14) 

where the two new normalized vectors are defined by 

| in) eSpan({|i(r))}), (15) 



|out)€Ker({|i(r))}). (16) 

The phase of these orbitals can be chosen so that the coefficients a and j3 are real and nonnegative, and 
a 2 + (3 2 = 1. We also rewrite the Slater determinant in terms of a new orthogonal basis of orbitals: 

niI 1 at(r)|0)=a| n n M at|0). (17) 

Here we introduce new orbitals | mu), fj, = 2, 3, N, such that Span({| i(r))}) = Span(| in), {| /i)}), that 
is, the space of filled states remains the same. Eq. (flTt can be shown by the following steps. Let U be a 
fermionic linear-optics transformation such that for i = 1, . . . N 

D 

tfajtft = Y, Vimcb, (18) 

m=l 

where the D x D unitary matrix V = V © W is a block-diagonal matrix with V G SU(iV). Furthermore V 
is such that | l(r)) is rotated to | in) and | i(r)), i > 1, is rotated to | /x), \i = 2, 3, N. Thus, as required, 
the rotation V does not change the space of filled orbitals. By inserting W U = I and using U\ 0) = | 0) we 
can see that 

C/n^ 1 a|(r)|0)=aJ n^lO), (19) 
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while we can also write 

N N N 

U Uf =1 a\ = Vu4 £ V 2j a] £ V 3k a\ ... (20) 

i=l j=l k=l 

There are N modes, and N sums in this expression. In order for a term to be nonzero, each mode must 
appear once and only once (otherwise a mode would be repeated, and (a)) 2 = 0). Thus we can write using 
the anti-commutation relations that 

U Ug ia \ = sign(7r)V 1)7r(1) V 2)7r(2) V 3 , 7r( 3 ) niI 1 at = det (V) D^a} = ug.a]. (21) 

Thus, Eq. (fT7b is established. 
Now we can calculate: 

P! (11^4(^1 0) = ata K (nf =1 aJ(r)|0) 

= 4(«Oin + / 3a out) a in II ^ a il ) =aa « n M4l°)- (22) 

Here we have used the fact that a ou t| 0) = and aj n a| n l 0) = | 0). We note that the final state is indeed just 
another Slater determinant; it is unnormalized, but the coefficient a just reflects the fact that the probability 
of this outcome is a 2 . Note that a is easy to calculate, as it is just the magnitude of the projection of a vector 
in the single-particle space of dimension D. The outcome of the other projector just takes a little bit more 
to evaluate: 

PoOl^K^IO) = a K at(n^ai(r)|0) 

= a K (aa| n + /?aJ )Ut )aJ n n At aJ t |0) = /3a K aJ )Ut a| n n At a];|0). (23) 

Here we have used («| n ) 2 = 0. To go further, we have to introduce another orbital | k -1 ) perpendicular to 
| k) in the space spanned by | in) and | out) ; specifically, 

| K x ) = j3\ in) - a | out). (24) 

Then using the relation 

a out a in = a « a !^> < 25 ) 

we finish the derivation: 

P (ldaJ(r)| 0) = /Ka^atn^at | 0) = ^al^TL^ | 0) = /^n^aj | 0). (26) 

So all operations keep the state vector in the Slater-determinant form. Indeed, we note that except 
for the change of normalization, the action of Po,i on the state is identical to that of some single-particle 
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Hamiltonian. This gives a new view on the classical simulatability of this suite of operations. Appendix A 
gives some correspondence between these calculations and ones that appear in the theory of electron energy 
bands in crystals. 

2 Two-mode measurements 

We now consider a different scenario for quantum computation, one in which one can perform a two- 
mode measurement. This can for example be a "charge measurement" that determines how many electrons 
are in a particular spatial orbital irrespective of their spin state. This means that we imagine a two-mode 
measurement, in which the two modes are identical spatially and differ only in their spin quantum number 
(so we might indicate these two orbitals | m j) and | m |), for some spatial orbital m). In the present 
analysis we will not wish to make the distinction between spin and orbital labels, so that we will just 
consider a number measurement in two orthogonal modes labelled | n) and | A). (This generalization means 
that our analysis is applicable to problems involving spin-orbit interaction, where the distinction between 
spin and space is not applicable.) 

We will assume that this measurement is "nondestructive", a feature of recent charge measurements in 
quantum dot structures [7] (but, see the remarks in the Discussion section). Then, similar to a single-mode 
measurement, we must write down the 'operation elements' (see [8 1) of the measurement which we take to 
be projectors (see the Discussion for potential modifications). 

For the "0" outcome (both modes unoccupied), this is simple, since it is just given by the product of the 
two "0" projectors for the two modes separately: 

Po = PokPoa = a K ala x a[. (27) 

The "2" outcome projector (both orbitals occupied) is also simply the product of the two one-orbital projec- 
tors: 

P2 = Pi^Pia = a{a K a[a x . (28) 
The "1" projector can be written as a sum of two projector products: 

Pi = a K a) K a\a\ + a[a K a x a[. (29) 

The important point for the upcoming analysis is 

Lemma 1 Pi cannot be expressed as a single term; that is, it is not possible to write P\ as P\ = /i/2/3---/m» 
where M is an arbitrary integer, and each fi is either a creation or an annihilation operator for some arbi- 
trary mode. 
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Proof: We have to study two cases: 

1) /m is a creation operator. We consider a one-particle basis consisting of orbitals | k) and | A), and 
D — 2 orbitals | u) orthogonal to | k) and | A). Then the orbital 4>m that Jm creates can be written 

\<f> M )=a\K)+0\X) + ^2c v \v). (30) 

Consider the (unnormalized) state 

|$) = (aal + pa{)TL v al\0). (31) 
This state has one electron in the space Span(| k), | A)), so it is an eigenstate of Pi with eigenvalue 1: 

Pi| $) = |$). (32) 
But Jm annihilates | <£) (recall that (at) 2 = for any x): 

f M \ $> = (aal + (3a{ + ^ cv4)("4 + M) n ^4l 0) 
= [(aat + Pa\) 2 U u ai - ^ c v ,{aa{ + Pa\)al,Tl„ai}\ 0) = 0, (33) 

u' 

so the single-term expression /i/2/3---/a/ cannot equal Pi in this case. 

2) /m is an annihilation operator. We consider the same orbital expansion as in Eq. d30b . and we apply 
/m to the state 

|*) = (a*o« + ^o A )44l > ( 34 ) 
This is again an eigenstate of Pi, Pi| 1 f) = | ^). However, 

f M \ V) = {a*a K + (3*a x + c* u ,a u/ )(a*a K + /3*a A )44l °> 

v' 

= [(a*a K + (i*a x fa) R a\ - ^ c*,(a*a K + /3*a A )a^4 a ^']l 0) = (35) 
(since (a x ) 2 = and a x \ 0) = 0). So, in this case too the single-term expression cannot match Pi. □ 

So, the fact that Pi cannot be written as a single term opens the possibility that this two-mode mea- 
surement can lead to more complex quantum time evolution, and thus has the possibility of implementing 
quantum computation. In particular, when the "1" outcome is obtained, the fact that the minimal expression 
for Pi contains two terms means that the Slater number (recall above, see [4|)) could double after every 
such measurement; so, if there are M such "1" outcomes, then the state may have an exponentially large 
(2 M ) Slater number, a state for which expectation values are likely to be very hard to calculate classically. 
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We do not know how to prove that the Slater number will in fact be as large as 2 M although we will 
now prove that generically, when Pi is applied to a single Slater determinant, the result has Slater num- 
ber two. Nevertheless, the expectation that the Slater number becomes high, and the evolution becomes 
difficult to simulate, is vindicated by the discovery of Beenakker et al. [ 1 ] that quantum computation is 
implementable by linear fermion optics and the two-mode measurement! We will return to more discussion 
of this measurement after Sec. [3] 

3 Slater number generically goes from one to two under Pi 

We now show that if we apply the two-mode projector Pi to a single Slater determinant for N > 2 electrons, 

| *) = Pin^aJl 0) = {a K ala\a x + (J K a K a x a\)B.? =l a\\ 0), (36) 

then the resulting state | \F) generically has Slater number two. Note that we can always choose a basis 
such that the initial state has the standard form shown. Furthermore, without loss of generality, we can 
parameterize the orthogonal orbitals | k) and | A) as 

\k) = cos#| 1) +sin#|iV + 1), (37) 
| A) = cos0(-sin0| 1) +cos0\N+ 1)) + sin 0(cos £| 2) + sin£| iV + 2}). (38) 

We can simplify the problem considerably by using a theorem of K. Eckert et al. @, that the Slater number 
cannot be increased by applying an annihilation operator to a state. Since | k) and | A) do not involve orbitals 
| 3), | 4), ... | N), we can annihilate all of these and be left with a two-electron state: 

| *') = Uf =3 ai\ = {a K aia[a x + ala K a x a\)a\a\\ 0). (39) 

Using the methods of Sec. ^ we can convert each term of this expression into one involving just two 
creation operators. After a lengthy calculation (using Mathematica) we find 

l^)= w ija \a)\Q). (40) 

i,j=l,2,N+l,N+2 

Where the 4x4 antisymmetric matrix w is 



^ D fa nnsfl fr~i sin f) cos§sin2</> cos 8 sin 2 <f> sin 2£ sin 6 sin 2 <j> sin 2£ N 

' S 2/ s 2/c 



fs cos 9 - f c sin B 

-f c cos 6 - f s sin 6* 





v • 



sin 2<f> sin § 

2fc 

cos sin 2 <f> sin 2g sin 9 sin 2 <f> sin2g 

Vc 2fi 




(41) 
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(we don't show the lower triangle of this antisymmetric matrix), with 



fc = \/ cos 2 4> + cos 2 £ sin 4>, 



fs = y cos 2 4> + sin 2 £ sin 2 <j>. (42) 

As discussed by [4|, a basis transformation can be made to bring an antisymmetric matrix to a canonical 
form, consisting of a direct sum of 2 x 2 antisymmetric blocks. The number of nonzero blocks is the Slater 
number. Obviously, for a 4 x 4 matrix the Slater number is two iff both blocks are nonzero, and iff the 
determinant of the matrix is nonzero. For an antisymmetric matrix it is more convenient to evaluate the 
Pfaffian, which is the square root of the determinant. For w we find that the Pfaffian is 

Pf(w) = (43) 

So, we see that generically, Pi does indeed increase the Slater number from one to two. 



4 A no-go theorem 

We now explore further the power of nondestructive two-mode measurements. As noted above, Beenakker 
et al. Q have shown that the two-mode electric charge measurement above, in conjunction with linear 
fermion optics, permits the efficient implementation of quantum computation. It was noted, however, that 
like most of the linear photon optics schemes proposed to date (cf. [9]), this implementation using the three- 
outcome charge measurement is non-deterministic, i.e., there is some finite chance that the computation 
fails (in the case of an unlucky combination of measurement outcomes), although the overall probability 
of failure can be made acceptably low by a suitable implementation strategy. We now argue that this small 
probability of failure is intrinsic to this implementation: 

Theorem 1 If there exists an efficient implementation of the unitary evolution of a quantum circuit using 
linear fermion optics ( including single-mode measurements) and the three-outcome, two-mode measurement 
of Sec. ^\that is exact, i.e., has zero probability of failure, then this unitary evolution has an efficient classical 
simulation. 

Proof: Suppose the exact implementation exists. The efficient classical simulation of this unitary evolution 
proceeds as follows: We begin with the standard, single Slater determinant state of Eq. Q. We compute 
the effect of the first stage of fermionic linear optics on this state as in Eq. @. Then, we consider the 
first 0/1/2 charge measurement. We can calculate whether the probability for outcome "1" is 100% or not 
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by computing the action of Pi on the state (a simple calculation for a single Slater determinant). If the 
probability of "1" is 100%, if the state is an eigenstate of Pi with eigenvalue one, then the measurement 
has no effect on the state, and we proceed on with the next stage of computation. If the probability of "1" is 
not 100%, then at least one of the outcomes "0" or "2" has nonzero probability. Since the implementation 
of the quantum gates is supposed to be exact (i.e. works for every measurement outcome) we are free to 
choose any outcome that occurs with non-zero probability. So we choose or 2 (making sure the choice 
has nonzero probability) and then note that the state after application of Po or P2 is still a single Slater 
determinant. By proceeding thus, the classical simulation at all stages need only keep track of a single 
Slater determinant, which is efficiently doable. □ 

Remark: In this proof we have assumed exact classical real-number computation. The proof can be 
relaxed to treat the case of finite precision classical computations; in that case the quantum computation is 
simulated approximately, but always with high precision. 

Thus if we believe (which we do) that there does not exist an efficient classical simulation of the uni- 
tary evolution of polynomially-sized quantum circuits 3 , it follows by this Theorem that there will be no 
exact implementation of quantum circuits using fermionic linear optics and the two-mode three-outcome 
measurement. 

We can modify the two-mode measurement such that some of the outcomes are not distinguished and see 
what happens. For example, we can consider a two-outcome measurement with projectors Po,i = Po + Pi 
and P2, which only distinguishes whether the two modes are completely filled or not. All such grouped 
measurements can function in the nondeterministic implementation of quantum computation of [ 1 1. But 

Corollary 1 Theorem[l\still holds if the three-outcome measurement is replaced by the two-outcome mea- 
surement Po,i/^2» or Pi^/Po- 

Proof: For both measurements there is an outcome (P2 in the first case, Po in the second) for which the 
simulated state remains a single Slater determinant. The rest of the proof then applies. □ 

However, for one measurement (the "parity" measurement), this argument does not apply: 
The no-go theorem does not apply to the parity measurement Pq^/Pi- 

It would apply if one of the projectors could be written as a single term. We have already demonstrated 

that Pi cannot be written as a single term. This is easy to show for the projector Po,2 = Po + P2 as 

well, by similar arguments: If /m (see Lemma 1) is a creation operator, it annihilates the I?-electron Slater 
3 Note this does not include the final single qubit measurements. 
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determinant, which is not annihilated by Po,2; if fu is a destruction operator, it annihilates the vacuum | 0), 
which is not annihilated by Pq^. 

It was this observation that led Beenakker et al. to investigate alternative implementations of quantum 
circuits using linear fermion optics using the two-mode parity measurement; and, indeed, an exact simu- 
lation, which is in some sense much more efficient than any of the known non-deterministic simulations, 
exists! 

5 Discussion 

We have presented an alternative description of the fermionic linear optics computation. It is likely that 
the extension to "fermion-parity preserving" quadratic Hamiltonians which was treated in Ref. 0, can be 
analyzed similarly using Slater determinants. 

We want to close with a few words of caution about the applicability of our results. We have indicated 
that the two-mode measurement that enables quantum computation is "nondestructive" and uses projective 
measurement 'elements'. What happens if we relax these conditions? 

If the measurement is destructive, it means that the modes | k) and | A) are no longer available for further 
processing. The "tracing out" of these two modes that this throwing away implies is implemented in second 
quantization in the following way: the density matrix of the system, after the application of the two-mode 
projectors discussed above, is changed by the application of two trace-preserving completely positive maps, 
T K and T\. The trace-over-^ map 7^ is given by 4 



This map leaves the one-mode measurements unchanged; but the two-mode measurements are changed in 
a very important way. In particular, for any of the two-mode measurements discussed above, the tracing 
out leaves the system in a density matrix that is a mixture of single Slater determinants. The evolution of 
such states is simple (that is, efficiently simulatable on a classical computer), so destructive measurements 
give none of the quantum computational power of nondestructive ones. This has been anticipated in earlier 
studies of quantum measurements for quantum computation [ lOl fTTl . 

Another important modification of the measurement is the following. Instead of the two-mode parity 
measurement with measurement elements 1 8 1 Po,2 and Pi, suppose we have a two-mode parity measurement 
4 Actually, the simpler choice of Kraus operators A\ — aj, Ai = has the equivalent effect. This corresponds to going to a 
hole representation for the mode 



2 




(44) 
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with measurement elements Ui n tPo,2, (i-e. a projector), and Pi, where j7» n t is a charge-preserving 
unitary interaction. The probabilities of outcome of these two measurements are the same but the state after 
the outcome 0/2 has occurred will have undergone an additional unitary transformation Ui n t in the second 
type of measurement. 

The status of the no-go theorems is the same for these two measurements, but the construction in Ref. 
[ 1 1 only applies to the first. 

This could be important, as there may be situations where U{ n t is nontrivial. In particular, the Beenakker 
construction Q is isomorphic to one in which the qubit is coded by one electron in a double quantum dot, 
with occupation of the orbital in the left dot representing | 0) and right-dot orbital occupied being | 1) [ 12]. 
Then, as is suggested in [1], the required charge parity measurement might be accomplished by placing a 
single-electron transistor between the right dot of one qubit and the left dot of the adjacent qubit, so that it 
is sensitive to the charge on both (and can be tuned so that it reads one value of current for both dots empty 
or occupied, and another level otherwise). However, an analysis of this setup 1131 might reveal that Ui n t is 
nontrivial in this case due to an effective interaction which the measurement sets up between the electrons 
in the two qubits. More analysis would be worthwhile in this case; and if, in this or in similar situations, 
Ui n t turns out to be different from the identity, it would be worthwhile to work out an implementation of 
quantum gates for this case. 

Acknowledgments 

We are grateful for the support of the National Security Agency and the Advanced Research and Develop- 
ment Activity through contract DAAD19-01-C-0056. 

Appendix A: Application to band theory 

The analysis of Sec. [2can be applied to simple model problems in electron band theory. Not surprisingly, 
given the long history of band theory, see e.g. fl4l . some of the objects obtained above have special names, 
and special significance, in this setting. 

Suppose we consider a non-interacting Hamiltonian for electrons on a one-dimensional lattice. If the 
Hamiltonian only contains nearest-neighbor hopping terms, iaja, + i, and t < 0, then the ground state of the 
system is the Slater determinant 

nfi< feF 4lo>- (45) 
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The orbital a\ | 0) = | k) is the plane wave 

\k) = ^=Y. elkx \ x )- < 46 > 

* X 

Here D is the number of lattice sites (we assume periodic boundary conditions), a£| 0) = | x) is the orbital 
centered at site x, (x = 0, 1, 2, D — 1), and k assumes the values k = 2irn/D (k lives in the reciprocal 
space of the crystal), with integer — (D — l)/2 < n < (D — l)/2 (The electrons are, in this example, 
spinless). The filled states have |/c| < 27r(iV — 1)/ (2D) = /cp where the Fermi-wavenumber fcp is 27r(iV — 
1) / (2D) ~ ttu for iV >> 1. We assume that iV is odd. Here v = ^ is the fiiiing of the band (the number 
of electrons per orbital | x)). Note that the empty states are those with kp < \k\ < ir, and k = tt, but not 
k = —it, corresponding to the rule that A;s differing by a reciprocal lattice vector, in particular those lying 
on the boundary of the first Brillouin zone, should not be counted twice. 

Suppose that a measurement is done that reveals that an electron is present at the origin. What is the 
new Slater determinant describing the state? That is, we are specializing the development in the text to the 
case 



k> = |z = 0> = J=£|fc>, (47) 



k 



l in > = 7^ E l fc ) = l^o), (48) 

\k\<k F 

|out) = - 1 V \k), (49) 

\k\>k F 

\ K ) = y/D\m) + y/\-v\ out), (50) 

\ W ») = -7Tt E etkS \ k ^ « = 0,1,...,JV-1. (51) 

VN \k\<k F 

Here we have introduced the orbitals | W s ), which are obtained by a Fourier transform over the filled states 
| k). They are somewhat localized on the lattice, but not perfectly, since they only include the plane waves 
up to a certain wavelength. The wave function of Wq is 

W (x) = (x\W ) = ^pl. (52) 

The other orbitals | W s ), s / 0, are basically displaced versions of Wq: 

W a (x) = Wo(x-s/v). (53) 

Note, however, that an analytic continuation of x to the reals is understood here, since for general v the W s 
wave functions are generally not centered on lattice sites. 
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The orbitals | W s ) rather resemble the Wannier functions of band theory, in that they are approximately 
localized states built out of band orbitals. They are different, though, in that Wannier functions are generally 
defined for full bands only, i.e., only for v = 1. 

So, again, how is the fermi sea modified if the electron number is measured at the origin? With proba- 
bility v the answer will be "1", and then the new Fermi sea has the Wo orbital replaced by the completely 
localized orbital | k) = | 0) , and all the rest unperturbed: 

4n^ <4j0). (54) 

One can say that one orbitals' worth of electrons has been pulled out from 0(1 jv) lattice sites around the 
origin and concentrated at x = 0. The hole that is left in the fermi sea by this process is what is known as 
the exchange hole in electron physics [ 15 ]. 

With probability 1 — v the measurement gives answer "0"; then the Wq orbital replaced by | « ): 

o^n^okjO). (55) 
This modified orbital can be rewritten in an informative way: 

| = VT^ | in) - 07 1 out) = -Jr^—l 0) + \ T^—\ Wo). (56) 

VI — v VI — v 

We find that the wave function for this orbital (x {k- 1 } is zero for x = 0, as expected, and also has an ex- 
change hole, but with a reversed sign compared with the other measurement outcome, and with a magnitude 
that depends on v. If v is near one, the perturbation of the fermi sea is very strong, but also this outcome 
"0" occurs with vanishingly small probability. 
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